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In wireless communications, spatial diversity techniques, such as space-time block code (STBC) and 
single-input multiple-output (SIMO), are employed to strengthen the robustness of the transmitted signal 
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t-^J- ' against channel fading. This paper studies the performance of protograph-based low-density parity-check 
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(LDPC) codes with receive antenna diversity. We first propose a modified version of the protograph 
extrinsic information transfer (PEXIT) algorithm and use it for deriving the threshold of the protograph 
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f~>) 1 codes in a single-input multiple-output (SIMO) system. We then calculate the decoding threshold and 

simulate the bit error rate (BER) of two protograph codes (accumulate-repeat-by-3-accumulate (AR3A) 
code and accumulate-repeat-by-4-jagged-accumulate (AR4JA) code), a regular (3, 6) LDPC code and 
two optimized irregular LDPC codes. The results reveal that the irregular codes achieve the best error 
performance in the low signal-to-noise-ratio (SNR) region and the AR3A code outperforms all other 
codes in the high-SNR region. Utilizing the theoretical analyses and the simulated results, we further 
discuss the effect of the diversity order on the performance of the protograph codes. Accordingly, the 
AR3A code stands out as a good candidate for wireless communication systems with multiple receive 
antennas. 
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Channel state information (CSI), extrinsic information transfer (EXIT) algorithm, protograph-based 
LDPC code, receive diversity, single-input multiple-output (SIMO) system. 
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I. Introduction 

In wireless communications, fading is a major factor that deteriorates the quality of signal transmission. 
Many methods have been proposed to mitigate the effect of fading. Of particular interest are the multi- 
antenna technologies which can provide high diversity gain and spatial multiplexing gain (TJ. Recently, 
a wealth of research has investigated the interplay between forward error correction (FEC) and spatial 
diversity. As a type of superior FEC code, low-density parity-check (LDPC) codes are known to perform 
near the Shannon limit over the additive white Gaussian noise (AWGN) channel [2]. However, LDPC 
codes that perform well in AWGN channels may not do so in a fading environment [3]. To overcome this 
weakness, LDPC codes have been studied in fading environments and capacity-approaching LDPC codes 
have also been designed for single-input multiple-output (SIMO) channels and multiple-input multiple- 
output (MIMO) channels using density evolution S, Q, O and extrinsic information transfer (EXIT) 
function (3,(11, respectively. 

Many of the capacity-approaching LDPC codes, however, are irregular and hence suffer from two main 
drawbacks — high error floor and nonlinear encoding. Recently, some research has successfully reduced 
the error floor of short-block-length LDPC codes which can then attain outstanding performance down to 
an error rate of 10~ 5 O- Meanwhile, the quasi-cyclic LDPC codes that permit linear encoding have been 
proposed for MIMO channels [ 1 ] . In addition, a novel class of LDPC code, namely multi-edge type 
(MET) LDPC code, has been introduced (TTI . As one subclass of the MET-LDPC code, the protograph- 
based LDPC code has emerged as a promising FEC scheme due to its excellent error performance 
and low complexity (121 . Two types of protograph codes, namely accumulate-repeat-accumulate (ARA) 
code and accumulate-repeat-by-4-jagged-accumulate (AR4JA) code, which can realize linear encoding 
and decoding have been proposed by Jet Propulsion Laboratory (JPL) (TBI , (141 . (T51 , (161 . In (171 , a 
protograph EXIT (PEXIT) algorithm has been introduced to predict the threshold of protograph codes 
over the AWGN channel. While the protograph codes have further been studied under Rayleigh fading 
channels (181 . to the best of our knowledge, little is known about the analytical performance for protograph 
codes under fading conditions and antenna diversity has never been considered. 

In this paper, we aim to investigate the performance of the protograph codes over a SIMO Rayleigh 
fading channel. To do so, we propose a modified PEXIT algorithm for analyzing the protograph LDPC 
code over a fading environment. We then analyze the decoding thresholds of the accumulate -repeat-by- 
3-accumulate (AR3A) code, the AR4JA code, the regular (3, 6) LDPC code and two optimized irregular 
LDPC codes f4] and use the thresholds to predict the error performance of the codes. The results show 
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that the irregular LDPC codes outperform the other codes in the low signal-to-noise-ratio (SNR) region. 
However, in the high-SNR region, the AR3A code possesses the best error performance. Besides, we 
also study the performance of the AR3A code and the AR4JA code with different diversity orders based 
on (i) the mean and variance of the log-likelihood-ratio (LLR) values, (ii) the PEXIT analysis, and (iii) 
the bit-error-rate (BER) simulations. We find that the additional gain becomes smaller as we increment 
the diversity order and hence increase the system complexity. 

We organize the remainder of this paper as follows. In Section [TTJ we present the system model over 
the SIMO fading channels and in Section [Till we describe our modified PEXIT algorithm for analyzing 
protograph codes in the fading environments. In Section JV] we analyze the decoding threshold and 
the initial LLR distribution of two conventional protograph codes. We show the simulation results in 
Section [V] and finally we give the concluding remarks in Section [VT] 

II. System Model 

The system model being considered in this paper is illustrated in Fig. Q] Referring to this figure, the 
information bits (Info bits) are firstly encoded by the (punctured) protograph LDPC code. Then the binary 
coded bits v € {0, 1} are passed to a binary-phase-shift-keying (BPSK) modulator, the output of which is 
given by x = (— l) v G {+1,-1}. The modulated signal x is further sent through a SIMO fading channel 
with one transmit antenna and receive antennas. 

We denote h as a channel realization vector of size x 1, the entries of which are complex independent 
Gaussian random variables with zero-mean and variance 1/2, i.e., M(0, 1/2), per dimension. Moreover, h 
is assumed to be independent over time flU, 13, Q, (SJ. We define n as a x 1 complex AWGN vector 
with zero-mean and covariance matrix E(nnt) = (iVo/2)Ijv R = c^Nr where E(-) is the expectation 
operator, < denotes the transpose conjugate operator, and Iyv R represents the A^r x ATr identity matrix. 
Then, the Afo x 1 receive signal vector, denoted by r, is given by 

r = hx + n. (1) 

Note that up to now, the time index has been omitted for clarity. 

At the receiver, we assume that the Afo received signals in r are combined using the maximum- 
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Fig. 1. System model of the protograph codes over the SIMO Rayleigh fading channels. 



ratio-combining (MRC) method |[T9l or the equal-gain-combining (EGC) method GUlf* . Afterwards, the 
combined signals y are sent to the LDPC decoder for finding the valid codewords. We do not consider 
interleaving in our system^. Moreover, we assume that (i) the channels formed by different transmit-receive 
antenna pairs are independent and (ii) each receive antenna possesses perfect channel state information 
(CSI) which is changing sufficiently rapidly to satisfy ergodicity Q. The capacity of an ergodic SIMO 
channel is further given by HI 



C = E (\og 2 det (l NR + |^hh^ ) bits/s/Hz, 



(2) 



where E s is the average energy per transmitted symbol and det(-) is the determinant operator. This 
capacity can be evaluated by using Monte Carlo simulation, i.e., by generating a large number of 
independent channel realizations and computing their average capacity value. In practical applications, 
we should choose a modulation technique and a code rate R such that the product of the modulation 
order M c and the code rate equals the channel capacity, i.e., M C R = C. For BPSK modulation, M c = 1 
and hence we have R = C. Moreover, the definition of E^/No (Ei is the average energy per information 
bit) used in this paper is the same as that used in Q, i.e., 

E b (N R /R)E S 



N 



N 



(3) 



'Minimum mean square error (MMSE) combiner is particularly useful in mitigating interference in the frequency selective 
channels. Thus, a MMSE combiner is suitable for the cases in which interference exists 1211 . In this paper, however, interference 
does not exist in the channel. Consequently, we apply the simpler MRC and EGC 1191 , 1201 , which have also been used to 
process the received signals over interference-free Rayleigh fading channels incorporated with multiple antennas and LDPC 
codes (4), (5). 

2 Note that in a practical environment, interleaving may be required to ensure that the channel gains for the code bits in the 
same codeword are independent. 
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III. Modified PEXIT Algorithm for Protograph Codes 

The conventional EXIT function (TJ has been proposed to better trace the convergence behavior of the 
iterative decoding schemes and to efficiently estimate the thresholds of different codes. However, it is 
known not to be applicable to the protograph codes Q. In ifTTl . a protograph EXIT (PEXIT) algorithm 
has been introduced to facilitate the analysis and design of protograph codes over the AWGN channel. 
In the following, we illustrate that the PEXIT algorithm, which works well on the AWGN channel, is no 
longer applicable to a SIMO Rayleigh fading channel. Then we modify the PEXIT algorithm for such a 
channel and use it for analyzing the protograph codes in our system. 



A. Assumption of the PEXIT Algorithm 

One important assumption of the proposed PEXIT algorithm in IfTTl is that the channel log-likelihood- 
ratio (LLR) messages should follow a symmetric Gaussian distribution. In the following, we briefly 
illustrate that this assumption cannot be maintained in the case of a SIMO Rayleigh fading channel and 
then we elaborate how to apply the PEXIT algorithm in such an environment. To simplify the analysis, 
we assume that the all-zero codeword is transmitted. 

By using j (j = 1,2,...) to indicate the coded bit number and k (k = 1, 2, ... , Nr) to indicate the 
receive antenna number, the signal of the j th coded bit at the fc th receive antenna can be written as 



rj[k] = hj[k]xj +rij[k]. 

The combiner output corresponding to the j th coded bit, denoted by yj, is then given by |fl9l , ll20l 



(4) 



h*[k]rj[k] for MRC 



Vj 



£■-. 

k=l 

Nr h*[k] 



(5) 



E 



< k=l 1 J 1 



Tj[k] for EGC 



where * denotes the complex conjugate, | • | represents the modulus operator and K[k]/\hj [k]\ is used to 
remove the phase ambiguity for coherent reception in EGC. Subsequently, the initial channel LLR value 
L c h j corresponding to the j th coded bit can be obtained using HI 



L c h, 




In 



Pr(xj = +l\yj,hj) 



Pr(xj 
for MRC 

for EGC 



k=1 ^ (6) 



0fc=l 



where Pr(-) is the probability function and hj = [hj[l],hj[2], • • • , hj[Nji]] T (the superscript "T" repre- 
sents the transpose operator). 

To evaluate the performance of the two combiners, we examine the distribution of L c h.j by exploiting 
Monte Carlo simulations. We use a rate-1/2 AR3A code with an information length per code block of 
1024. Moreover, we consider a SIMO Rayleigh fading channel with N R = 2 and E b /N = 2.6 dB. By 
sending Xj = +1 repeatedly while varying the channel fading vector hj from bit to bit, we evaluate 
the mean of the absolute value of L c hj- We observe that the MRC produces an average value of 3.628, 
i.e., IEmrc(| Aihj'l) = 3.628 whereas the EGC gives EegcCL^c/ijI) = 3-242. Moreover, both combiners 
produce channel LLR values with almost the same variancdj. Consequently, using MRC should provide 
a higher chance for successful decoding. For this reason, we will focus on MRC in the sequel. 

We denote the real part of L c h by L re j, i.e., L re j = ^[L^ ,•]. In Fig.[2l we further plot the probability 
density function (PDF) of the L re j values (denoted by f(L re j)) when MRC is used. We also define 
uq = EMRc(^re,j) and Cq = var(L re j) and plot the PDF of the Gaussian distribution J\f(uo, ctq) in the 
same figure for comparison. The curves in the figure indicate that the PDF of the L re j values does not 
agree with the PDF of M(uq, On), which further suggests that the L c hj values do not follow a symmetric 
complex Gaussian distribution.]. We conclude that the PEXIT algorithm in ifTTl is not applicable to this 
type of channel. In the following, we analyze the distribution of the L^j values when the channel 
realization is fixed. 

We consider a fixed channel realization, i.e., a fixed channel fading vector hj. We assume using the 
all-zero codeword (i.e., xj = +1) and we substitute (0]) into ©. Then, we can rewrite the expression for 

L>ch,j as 



2 Nr 

n k=l 

2 ^ 



— E (\ h M 2 + h*[k} nj [k} . (7) 



3 The variance of a complex variable z is given by var[z] = E[(z — E(z))(z — E(z))*]. 
4 The same observation is found when EGC is used. 
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Fig. 2. Probability density functions of the L re j values and N(uq,o-q) over the SIMO Rayleigh fading channel. 



As E ( nj[fe]n*[/c] J = <r„, we have 



\hM l + h*[k] nj [k] ~ CM ( \hj[kr , Mm 2 °l 



and hence 



where 



\ n k=l n k=l J ^ n 



a 



2 a J 



fc=l 



(8) 



(9) 



(10) 



is defined as the channel factor. In this case, we find that the channel LLR values follow a symmetric 
complex Gaussian distribution. 

Besides, assuming E s is normalized to 1, ([3]) can be simplified to E^/Nq = (Nr/R)/No. Thus, 
combining it with iVo = 2<r^ gives 



N R 



(ID 



2R(E b /N ) 

To summarize, we observe that the channel LLR values for a SIMO fading channel do not follow a 
symmetric Gaussian distribution and hence the PEXIT algorithm cannot be applied directly. However, 
for a fixed fading vector, these LLR values follow a symmetric complex Gaussian distribution. Using this 
property, we propose a modified PEXIT algorithm that can be adopted to analyze the protograph codes. 
Details of the algorithm are described as follows. 
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B. Modified PEXIT Algorithm for Protograph Codes 

We first define some symbols and terms. A protograph G = (V, C, E) consists of three sets V, C and E 
corresponding to the variable nodes, check nodes and edges, respectively lfl2l . In a protograph, each edge 
eij G E connects a variable node Vj G V to a check node q G C. Moreover, parallel edges are allowed. 
A large protograph (namely a derived graph) corresponding to the protograph code can be obtained by a 
"copy-and-permute" operation. Hence, codes with different block lengths can be generated by performing 
the "copy-and-permute" operations different number of times. A protograph with N variable nodes and 
M check nodes can be represented by a base matrix B of dimension M x N. The (i, j) th element of 
B, denoted by b^j, represents the number of edges connecting the variable node Vj to the check node 
Cj. We also define five types of mutual information (MI) as follows. 

1) lAv(hj) denotes the a priori MI between the input LLR value of Vj on each of the bij edges and 
the corresponding coded bit Vj. 

2) lAc(hj) denotes the a priori MI between the input LLR value of a on each of the b{j edges and 
the corresponding coded bit Vj. 

3) lEv(hj) denotes the extrinsic MI between the LLR value sent by Vj to Cj and the corresponding 



4) lEc{i,j) denotes the extrinsic MI between the LLR value sent by q to vj and the corresponding 



5) Iapp(j) denotes the a posteriori MI between the a posteriori LLR value of Vj and the corresponding 



In addition, during each iteration in the PEXIT algorithm, we have lAcihj) = lEv(i,j) and lAv(hj) = 
lEc(hj)- We also denote the maximum number of iterations in the algorithm by T max . Besides, we define 
two new terms called indicator function and punctured label. 

Definition 1. We define the indicator function <£(•) of an element b^j in the base matrix B as 



Hence, $(6ij) indicates whether Vj is connected to q or not. 

Definition 2. We define the punctured label Pj of a variable node Vj as if Vj is punctured, and 1 
otherwise. 



coded bit Vj. 



coded bit Vj. 



coded bit Vj. 




(12) 



Moreover, the MI between a coded bit and its corresponding LLR value L c h ~ j\f(^,a 2 h ) is denoted 
by J(a c h) and is expressed as Q 

,« exp(-ii^) 
J(cr ch ) = l- V - ' log 2 [1 + exp(-e)] de (13) 

v 27rcJ ^ 

The inverse function of ([T3l is further given by Q 

[ l\x 2 + 12% + l?,yfx if < x < 0.3646, 
J (x) = < (14) 
I 74 ln[75(l — x)] + 76^ otherwise, 

where 71 = 1.09542, 72 = 0.214217, 73 = 2.33737, 74 = -0.706692, 75 = 0.386013 and 76 = 1.75017. 

Then, for a rate-i? protograph with N variable nodes and M check nodes, the proposed modified 

PEXIT algorithm over a SIMO Rayleigh fading channel can be described as follows. 

1) For a given SIMO channel realization h = [h[l], h[2], ■ ■ ■ ,h[Nf(]] T , we can calculate the corre- 
sponding channel factor a using (TTOb . i.e., a = J2k=i IM^]| 2 - Suppose we are given the number 
of blocks of channel factors (denoted by Q) and the maximum number of iterations (T max ). We 
generate a matrix a = (a q j) = (%2k^i l^gjMI 2 ) °f dimension Q x N to represent the Q blocks of 
channel factors, i.e., each row in a. represents a group of channel factors for the N variable nodes 
in the protograph. We also select an initial E^/Nq (in dB) which should be sufficiently small. 

2) For i = 1, 2, . . . , M and j = 1, 2, . . . , N, we set the initial lAvihj) to 0. We also reset the iteration 
number t to 0. Considering the punctured label and substituting (|TTb into (©, for the channel factor 
&q,j (j = 1,2, ... ,N and q = 1,2,..., Q), the corresponding variance of the initial LLR value 

'ch,q,j) 



(denoted by a 2 ch ) is given by 



2 _ ^" 3^1,3 



4P,Q! A &RPjOL a j „ _ (Efc/jVp) 



^10 • . (15) 



cK ™ ol N R 
3) If t = T max , set E b /N = E b /N + 0.001 dB and go to Step El otherwise, for i = 1,2,..., M; 
j = 1, 2, . . . , N and q = 1,2,...,Q, we calculate output extrinsic MI sent by Vj to q for the q th 
fading block using ifTTl 

/ 



lEv,q(hj) = ®{h,j)J 



(16) 



4) For % = 1, 2, . . . , M and j = 1,2, . . . ,N, we obtain the expected value of lEv,q(hj) using 

Q 



WEv,q{i,j)] = ^Y,lEv,q{i,j). (17) 
9=1 
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Then, the a priori MI between the input LLR of Cj on each of the 6jj edges and the corresponding 
coded bit is evaluated using 

I Ac {i,j)=WEv, q {hj)]. (18) 

5) For i = 1,2, ... ,M and j = 1,2, . . . , N, we compute the output extrinsic MI sent by Cj to vj 
using H71 

/ 



/ £c (i,i) = *(M i- J 



Vi 



E^fJ-Hl-^M))] 2 +(6y-l)[J- 1 (l-/Ac(«,j))] S 

(19) 



Then, we get the a priori MI between the input LLR of Uj on each of the bi j edges and the 
corresponding coded bit using 

lAv(i,j) = Ie&J)- (20) 
6) For j = 1, 2, . . . , N and q = 1, 2, . . . , Q, we compute the a posteriori MI of Uj using ifTTl 



Wi) = J U (E M^OM*, j))] 2 ) + <^ lW - j • (2D 

Then, for every j = 1, 2, . . . , N, we can evaluate the expected value of I app ,q(j) using 

1 Q 

E^opp.gCj')] = qY1 ^PpM)- (22) 
g=l 

7) If the expected MI values E[I apPt q(j)] = 1 for all j = 1,2, ... ,N, the E b /N value will be the 
EXIT threshold that allows all variable nodes to be decoded correctly and the iterative process is 
stopped; otherwise, we increase t by 1 and go to Step [3] to continue the iterative process. 

Note also the following. 

• Each of the output extrinsic MI sent by vj to q depends on both the initial channel LLR values 
&ch qj < n3T ) AND the type of code used. Thus, for a fixed we can compute Q different values 
for lEv,q(i,j) CED- Consequently, an average quantity E{lE v ,q{i, j)] (fT71 ) can be obtained based on 
the Q different values of lEv,q(hj)- 

• The average quantity K[lE V ,q{h j)] becomes the a priori MI I^ c (i, j) between the input LLR of c, 
on each of the b^j edges and the corresponding coded bit ( fT8l ). 

• Since the output extrinsic MI sent by q to Vj is only related to the type of code used and is 
independent of the initial channel LLR values o~ 2 chq -, the output extrinsic MI sent by Cj to vj, 
i.e., lEc{i,j), can be computed using ( fT9l , which depends only on the a priori MI I Ac (i,j) (i = 



1 1 



1, 2, . . . , M and j = 1, 2, . . . , N) and is independent of a 2 ch ■. For a fixed (i, j), there is only one 
lEc(i,j) computed and thus no averaging is required. 

The quantity lEc{hj) then becomes the a priori MI lAv{hj) between the input LLR of Vj on each 
of the bij edges and the corresponding coded bit (l20b . 

Each of the a posteriori MI value I a pp,q(j) of Vj depends on both the initial channel LLR values 
a ch q j AND the type of code used. Thus, for a fixed j, we can compute Q different values for 
I app ,q(j) (I2TT) . Afterwards, an average quantity E[I apP!g (j)] (l22l can be obtained based on the Q 
different values of I app ,q(j)- 

To ensure the accuracy of the modified PEXIT algorithm, we should generate a sufficiently large 
number of blocks of channel factors, i.e., a large value for Q. In this paper, we use Q = 10 5 . 



IV. Analysis of Protograph Codes 

Protograph codes not only enable linear encoding and decoding to be implemented easily, but also have 
superior error performance over the AWGN channel lfT2l . As two typical LDPC codes constructed by 
protographs, the AR3A code and the AR4JA code possess excellent performance in the waterfall region 
and the error floor region, respectively, over the AWGN channel fi"3l . lfl4l . lfl5l . lfl6l . The corresponding 
base matrices of the AR3A code and the AR4JA code with a code rate R = (n + l)/(n + 2) are denoted 
by Ba3 and T$A4, respectively, where 



/ 



2n 



B 



A3 



\ 



121000 ••• 
021112 1 ••• 2 1 

\0 1 2 1 1 1 2 ■•• 1 2/ 

In 



(23) 



B 



Ai 



(24) 



120000 ••• 
031113 1 ••• 3 1 
\0122113 ••• 13/ 

We assume that the j th column of the matrix corresponds to the j th variable node and the i th row of 
the matrix corresponds to the i th check node. Note that the variable nodes corresponding to the second 
columns in (l23l and (l24l) are punctured. 

Using the modified PEXIT algorithm proposed in Sect. IIII-B I we firstly investigate the decoding 
thresholds of the AR3A and AR4JA codes with a code rate of R = 1/2 (i.e., when n = 0). For 
comparison, the regular (3, 6) LDPC code and the two optimized irregular LDPC codes (denoted as 
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TABLE I 

Decoding thresholds (E b /N )th (dB) and the capacity gaps A of the AR3A code, AR4JA code, regular 

(3, 6) CODE, AND THE OPTIMIZED IRREGULAR LDPC CODES (IRREGULAR CODE A FOR JVr = 2 AND IRREGULAR CODE B 
FOR Nr = 4) OVER THE SIMO RAYLEIGH FADING CHANNELS. THE PARAMETERS USED ARE n = AND R = 0.5. 



Nn 


Capacity 


AR3A code 


AR4JA code 


Regular code 


Irregular code in (4) 


(E b /N ) th 


A 


(E b /N ) th 


A 


(Et/N ) th 


A 


(E b /N \ h 


A 


2 


-0.514 


1.258 


1.772 


1.433 


1.947 


1.993 


2.507 


1.175 


1.689 


4 


-0.662 


0.871 


1.533 


1.011 


1.673 


1.535 


2.197 


0.823 


1.485 



irregular LDPC code A and irregular LDPC code B) in IH are used to gauge the performance. Moreover, 
the degree distribution pairs of the irregular codes are given as 

X A (x) = 0.270234x + 0.266315x 2 + 0.463451x 9 

(25) 

PA (*) = 0.566545x 6 + 0.433455x 7 
and 

A B (x) = 0.285637x + 0.285602x 2 + 0.428724x 8 

(26) 

Pb(x) = 0.998857x 6 + 0.001143x 7 
As shown in H, irregular code A and irregular code B are optimized for the diversity orders of Ar = 2 
and Ar = 4, respectively. For our system, the diversity order equals Nr. Table U shows the decoding 
thresholds and the capacity gaps of these four codes over SIMO Rayleigh fading channels with two 
different diversity orders. Results in this table indicate that the decoding thresholds of irregular code 
A and irregular code B are smallest over such channels with the diversity order 2 and 4, respectively. 
These small thresholds suggest that the irregular codes should possess relatively better performance in 
the low-SNR region Ifl4"1 . |fT51 . However, the irregular codes may be outperformed by other codes in the 
high-SNR region because of the error-floor issue. Also, the results demonstrate that the regular LDPC 
code possesses the highest threshold. Hence, the regular code is expected to be inferior in the low-SNR 
region. Between the AR3A code and the AR4JA code, it is further observed that the AR3A code has 
lower thresholds for both the cases of Ar = 2 and Ar = 4. 

We further investigate the performance of the AR3A code and the AR4JA code in such systems with 
different code rates and different diversity orders. The results are shown in Table ITT] Given a fixed code 
rate, we observe from the table that the decoding threshold of the AR3A code is smaller than that of 



5 The capacity gap is defined as the distance between the channel capacity and the decoding threshold. 
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TABLE II 

Decoding thresholds (Eb/N ) th (in dB) of the AR3A code and AR4IA code with different code rates R 

OVER THE SIMO RAYLEIGH FADING CHANNELS WITH DIVERSITY ORDERS JVr. = 1, 2, 3, AND 4. 



Code Rate 


AR3A code 


AR4JA code 


Nr = 1 


N R = 2 


N R = 3 


AT R = 4 


Nr = 1 


AT R = 2 


Nr = 3 


AT R = 4 


1/2 (n = 0) 


2.143 


1.258 


1.031 


0.871 


2.303 


1.433 


1.151 


1.011 


2/3 (n = 1) 


3.898 


2.540 


2.105 


1.903 


4.134 


2.752 


2.293 


2.084 


3/4 {n = 2) 


5.154 


3.402 


2.847 


2.609 


5.413 


3.624 


3.068 


2.783 


4/5 (n = 3) 


6.175 


4.137 


3.479 


3.156 


6.506 


4.322 


3.648 


3.315 


5/6 (n = 4) 


6.985 


4.698 


3.960 


3.616 


7.286 


4.871 


4.117 


3.746 


6/7 (n = 5) 


7.667 


5.129 


4.361 


3.969 


7.996 


5.294 


4.490 


4.103 


7/8 (n = 6) 


8.399 


5.559 


4.691 


4.292 


8.644 


5.691 


4.830 


4.406 



TABLE III 

The mean and variance of the initial LLR values output by the MRC combiner over the SIMO RAYLEIGH 

FADING CHANNELS. AN AR3A CODE WITH n = AND R = 1/2 IS USED. THE INFORMATION LENGTH PER CODE BLOCK IS 

1024 and E b /No = 2.6 dB. Diversity order N r = 1, 2, 3, and 4. 



Nr 


1 


2 


3 


4 


E(|i«*j|) 


3.615 


3.628 


3.634 


3.636 


va,r(L chJ ) 


20.251 


13.053 


11.669 


10.502 



the AR4JA code for all diversity orders under study. We also find that as the diversity Ar increases, the 
threshold decreases and hence the error performance of both codes should improve in the low-SNR region. 
The decrease in the threshold, however, is reduced as Ar is incremented. For example, we consider the 
AR3A code with a rate of R = 1/2. The threshold reduces from 2.143 to 1.258 as Ar increases from 
1 to 2; but further reduces to 1.031 and 0.871 only when Ar increases to 3 and 4, respectively. 

For the AR3A code, we also evaluate the mean and the variance of the initial LLR values output 
by the MRC combiner using Monte Carlo simulations. The results, which are generated with the same 
parameters used for creating Fig. [2l are given in Table [Hi] The results show that both (i) the rate that 
the mean of \L c hj\ increases and (ii) the rate that the variance of L c hj decreases are reduced as Ar 
is incremented. Similar simulations have further been performed for the AR4JA code and the same 
observations are obtained. These results imply that the rate of improvement of the error performance 
reduces when we increase the diversity order. 
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Eb/No(dB) Eb/No(dB) 



(a) (b) 
Fig. 3. BER and FER curves of the AR3A code, the AR4JA code, the regular LDPC code, and the irregular LDPC codes over 
the SIMO Rayleigh fading channels with the diversity order (a) — 2 and (b) — 4. 



V. Simulation Results 

In this section, we simulate the error performance of the AR3A code, the AR4JA code, the regular 
(3, 6) code, and the optimized irregular LDPC codes in flU over SIMO Rayleigh fading channels. We also 
discuss the influence of the diversity order on the error performance of the protograph codes. The code 
rate and the information length of each code block are R = 1/2 and 1024, respectively. We terminate the 
simulations after 500 bit errors are found at each E^/Nq. Also, the LDPC decoder performs a maximum 
of 100 BP iterations for each code block. 

A. Performance Comparison among Four Different Codes 

Figure [3] plots the bit error rate (BER) and frame error rate (FER) curves of the codes over the SIMO 
Rayleigh fading channels with diversity orders Ar = 2 and Ar = 4. It can be seen that the AR4JA code 
and the regular (3, 6) code are the two worst-performing codes for the two diversity orders. Moreover, 
referring to Fig. Ha), at a BER of 10~ 3 , the irregular code A has a gain of about 0.1 dB over the AR3A 
code, which remarkably outperforms the other two codes. Yet, the BER and the FER performance of the 
irregular LDPC codes has little improvement when the E^/Nq exceeds 2.4 dB upon which the error floor 
emerges. In the same figure, we observe that the AR3A code has excellent error performance for the 
range of E^/Nq under study. For instance, at E^/Nq = 2.8 dB, the AR3A code accomplishes a BER of 
4 x 10 -6 , while the irregular code, AR4JA code, and the regular code achieve BERs of 2 x 10 -5 , 10 -4 , 
and 8 x 10~ 4 , respectively. We also observe that at a BER of 10~ 5 , the AR3A code has a performance 
gain of 0.3 dB over the irregular LDPC code and the AR4JA code, which are superior to the regular 
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code. Moreover, a large gain can be expected at a lower BER. Similar conclusions can be drawn from 
Fig. [3fb), where the diversity order equals 4. 

In general, among the four types of codes, we conclude that the optimized irregular codes possess 
the best error performance in the low-SNR region while the AR3A code can provide excellent error 
performance in the high-SNR region over the SIMO Rayleigh fading channels. 

B. Discussion about the Diversity Order 

Figure |4] shows the BER results of the AR3A code and the AR4JA code for various diversity orders 
(Nr = 1, 2, 3, and 4). As can be seen from this figure, the AR3A code outperforms the AR4JA code 
by more than 0.3 dB at a BER of 2 x 10~ 5 for all diversity orders. We also observe that for a fixed 
BER, the required E^/Nq decreases as the diversity order increases. However, the rate of decrease is 
reduced with the diversity order. For example, we consider the AR4JA code at a BER of 2 x 10 -5 . The 
required E^/Nq decreases from 4.2 dB to 3.0 dB, 2.5 dB and 2.3 dB, respectively, as Nr from 1 to 
2, 3 and 4. Similar observations are found for the AR3A code. Furthermore, these observations agree 
well with the analytical results found in Section |IVl On the other hand, the system complexity increases 
with the diversity order. Accordingly, one should appropriately select the number of the receive antennas 
for a practical system so as to make a good balance between system performance and implementation 
complexity. 

10 \ 

10 2 

£io 3 

10 ' 



10 6 

1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 3.4 3.6 3.8 4 4.2 
Eb/No(dB) 

Fig. 4. The BER results of the AR3A code the AR4JA code over the SIMO Rayleigh fading channels with diversity orders 
N R = 1, 2, 3, and 4. 
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VI. Conclusions 

In this paper, we have studied the performance of protograph-based LDPC codes for SIMO systems 
under fading channel conditions. We have proposed a modified PEXIT algorithm for analyzing such 
systems equipped an MRC combiner. We have also used the proposed algorithm to evaluate the decoding 
threshold of the protograph codes and hence to analyze their error performance. Furthermore, we have 
compared the decoding thresholds and the distribution of the initial LLR values of the protograph codes 
among different diversity orders. We conclude that while the error performance of the protograph codes 
improves with the diversity order, the rate of improvement is reduced as the diversity order becomes larger. 
Thus, we have to strike a balance between code performance and system complexity when determining 
the diversity order to be implemented. Our simulation results have also shown that the AR3A code is able 
to provide a significant gain over the AR4JA code, the regular (3, 6) code, and the optimized irregular 
LDPC codes in the high-SNR region. In the future, we will strive to optimize the protograph codes in 
such systems. Furthermore, we will explore extending the modified PEXIT algorithm to other systems 
such as the MIMO fading systems. 
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